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Abstract: We consider the problem of optimal location of a Dirichlet region 
r) • in a d-dimensional domain Q subjected to a given right-hand side / in order to 

Q ■ minimize some given functional of the configuration. While in the literature the 

Dirichlet region is usually taken d — 1 dimensional, in this shape optimization 
problems, we consider two classes of control variables, namely the class of one di- 
mensional closed connected sets of finite one dimensional Hausdorff measure and 
the class of sets of points of finite cardinality, and we give a necessary condition 
of optimality. 
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1 Introduction 

We consider the problem of finding the optimal location of a Dirichlet region S 
in a d- dimensional domain Q associated the p-Laplacian 

'-diY{\Vu\P-^Vu) = f infi\S 
n = on afi U S, 

where the right hand side / is given as a nonnegative element of L^ {Q), being p' 
the conjugate exponent of p. The functional J-" we consider as a cost is defined 

by 

J-'(S) = / F(^x,u{x),\/u{x)) dx 
Jq 

where u is the unique solution of equation ([1]) and FifixIRxM"^— j-Misa 
Caratheodory function. The shape optimization problems we consider consists 
in the minimization of the functional J-" over two classes of admissible control 
variables S. The first class consists of all closed connected subsets T, oi Q whose 
one- dimensional Hausdorff measure 'H^(E) (sometime called the total length of S) 
is uniformly bounded by some constant L, while the second is the class of discrete 
subsets of fl whose 0-dimensional Hausdorff measure (i.e. their cardinality) does 
not exceed a given number n. 

The existence of an optimal configuration for the two optimization problems 
described above can be obtained by using a generalization of the Sverak's com- 
pactness result (see |Tl] for p = 2 and [3j for a general p). In this paper, we 
consider a penalized version of these shape optimization problems (see ([2]) and 
(12])) by adding to the cost functional J-" a Lagrange multiplier penalization of 
the form A'H^(S) and A'H°(S) respectively. This problem has been considered in 
[1] in the simplest case when p = 2, where the PDE ([1]) reduces to the classical 
Laplace equation. Moreover, it has been shown (see [5j) that, for the particular 
case where F{x,u,C,) = f{x)u, the limit problem, as p — )■ +oo, is the minimum 
problem for the average distance functional 



J^{^)= / dist{x,i:udn)f{x)dx, 

Jn 

where dist(x, E) denotes the distance between the point x and the set E. 

For the well posedness of the minimum problems ([2]) and ([3]) the constraint 
p > d — ds, where d^ is the dimension of S, has to be imposed in order to have 
Dirichlet regions of positive p-capacity. Let us mention that the regularity of the 
optimal sets (in the case of problem ([2])) is still an open problem, even in the 
simplest two dimensional setting with p = 2 and F = f{x)u. 



2 Setting of the problem and existence of mini- 
mizers 

We consider the following classes of control variables: 

A{fl) := {T.CQ : S closed connected, ^^(S) < +00} 

B{n) := {S C fi : S discrete, H\^) = #(S) < +cx)} 
and the shape optimization problems 

min{j^(S) + A-H^(S) : S G ^(fi)} (2) 

min {J^(S) + A-H°(S) : S G i3(fi)}. (3) 

The penalization terms A'H^(S) and A'H°(S) with A > replace the constraint on 
T-L^CE) and 'H*'(S) and prevent the minimizing sequences to spread over all the 
domain Q and hence getting a trivial solution. The existence of minimizers in 
the two shape optimization problems ([2]) and ([3]) is a consequence of the Sverak 
continuity-compactness result (see [llj for p = 2 and [3j for general p) and the 
Blaschke and Golab theorems. For the convenience of the reader let us give some 
details on the existence of an optimal shape. Let {S„}„ C A{Q) be a minimizing 
sequence in the optimization problem ([2]), then there exists a constant C > such 
that sup„{J-'(S„) + Ti^CEn)} < C. Since {S„}„ is a sequence a closed connected 
subsets of Q such that sup„?^^(E„) < C, by Blaschke theorem (compactness 
of the sequence {S„}„ in the Hausdorff topology) and by Golab theorem (lower 
semicontinuity of the "H^ with respect to the Hausdorff topology) , up to extracting 
a subsequence, {S„}„ converges in Hausdorff distance to some S G A{fl) and 
^^(S) < liminf„^oo^^(Sn)- For the lower semicontinuity of the energy part of 
the functional we need the Sverak continuity-compactness result which is stated 
in this terms: let {f2„}„ be a sequence of open and bounded sets contained in a 
fix bounded set D. If we assume that the number of the connected components 
oi D\fin is uniformly bounded by some number k, then {fln}n converges in the 
Hausdorff topology to some open and bounded set Q G D and the number of the 
connected components oi D \Q is less or equal to k. Moreover, if we denote by 
Un G Wf^'^i^Qn) the distributional solution of the p-Laplace equation 

-ApUn = f in fi„ 

u^ G W^o'''(^n) 

for some / in W~^'^ (D), then up to subsequence, «„ converges strongly in 
W^'P{D) {un are extended by zero outside i7„) to the function u which is the 
distributional solution of the equation 

-ApU = f in f2 
u G W^'^iQ). 



This result is interesting only in the case where p satisfies d — 1 < p < d because 
the case where p > d is trivial due to the fact that functions in W^'P{D) are 
continuous and the convergence of solutions follows easily. To apply this result to 
our problem, we choose Qn = -D\S„ and notice that {f2„}„ converges to fi = D\S 
in the Hausdorff topology where S is the limit of S„. From the hypothesis of the 
function F and the continuity with respect to the domains variation of solutions, 
the lower semicontinuity follows easily and also the existence of an optimal shape. 
The existence of an optimal set in the problem ([2]) is even easier due to the fact 
that p > d. 

Our goal is to derive first order necessary conditions of optimality, assuming 
that the solutions S of the minimum problems are regular as necessary. Notice 
that, since every set in S G A{Q) is countably "H^ rectifiable, it can be written 
as an union of countable C^ curves and a set of 7/^ measure zero. 

We distinguish three cases according to technical computations. We consider 
the following classes of control variables S: the class of closed connected subsets 
of R^, the class of points of W^ ioT d > 1, and the class of closed connected subsets 
of M'^ for (i > 2. In the two last cases, some extra difficulties occur because of the 
co-dimension of E which is greater than 1. For simplicity, we assume that Q has 
a smooth boundary and m = on dQ U S; we also assume as much as needed the 
regularity on the data. Before looking for the necessary conditions of optimality, 
we recall some definitions and results which will be helpful; we refer to [2] for the 
details. 

For a measure /i we denote for /i a.e. x by P^{x,-) : M'^ — )■ Tan(/i,x) the 
orthogonal projection of M'^ on Tan(yU, x) where Tan(/i, x) stands for the tangent 
space of /i at X that is the set of all tangent measures to n at the point x (see for 
instance [1]). 

Definition 1 The curvature of fi is defined as the vector valued distribution 

H^ := div(P^/i). 

In other words H^ is defined by 

{H^,X) = - [ div^ Xdfi VX e C^ (M^) , W^) , 



where div^X = E'=i(^m(VX^)),. 

We denote by A4bc the set of all positive and finite Borel regular measures 
of W^ whose curvature is a Borel regular measure with finite total mass. Since 
the curvature H^ of a measure /x G Mbc is not necessary absolutely continuous 
with respect to /x, by Radon-Nikodym theorem, we can write 

Hfj, = h{fi)fi + 9/i, 



where h{fi) G L)j(]R'^, M*^) is the density of iJ^ with respect to yU (also called the 
pointwise curvature) and dfi is the singular part of H^ with respect to /i (also 
called the boundary of fi). 

If /i = Ti^L.Ti, with S a C^ fc-manifold with boundary in M°', then by classical 
divergence theorem we have 

where h stands for the mean curvature vector of S and v the co-normal unit 
vector of dH. When the tangent space to yU is reduced to zero /i a.e., H^ is zero. 
This is for instance the case where /x is a finite sum of Dirac masses, or /i is 
concentrated on an a-dimensional Cantor subset C of [0, 1] with < a < 1 and 
?/"(C) G (0,+oo). 

Definition 2 Let T, be a countably T-i^ rectifiable set and fi = Ol-i^iSL be the asso- 
ciated rectifiable measure. A function h G L^(S,]R'^) is said to be the generalized 
mean curvature ofE if 

div^ Xdfi = - I X -hdfi WX e C^(R'^, R'^). 

In this case we denote the generalized mean curvature of S by Hj^. 

Theorem 3 Let {fir)r be a bounded sequence in M.bc weakly converging to some 
measure fi and assume that dimTan(/ir)/ir weakly converges to g^. Then the 
condition 

dimTan(/i,a;) < (^(a;) fx — a.e. (4) 

is necessary and sufficient to have 

P^^Hr -^ P^,^J'■ (5) 

In this case we have 

H,. - H,. (6) 

Proof: see [2j D 

Let n be an open subset of W^ and F : (] x R x M'^ -^^ [0, +oo] be a positive 
Caratheodory function. We assume F smooth and satisfying the growth condition 

F{x,u,z) < a{x) + \u\'P+ \z\P, 

where a is an L^{Q) function. We consider the functional 

J-'(S) := / F(x,'u(x), V-u(x)) (ix + Am(S) 

where m{Il) is either ?^^(E) if S is a closed connected one dimensional set or 
7/*'(S) if S is a discrete set of points and u the solution of the equation ([T]). We 
are interested in the necessary conditions of optimality satisfied by the minimizers 
of J-". From now on we assume optimal sets in the case of closed connected sets 
to be of class C^'° for some < a < 1 that is locally graph of C^'° functions. 



(7) 



3 Case of closed connected subsets in M' 

Let u be the weak solution of the state equation 

div (|Vm|p-2Vm) = / infi\S 
M = on afi U S, 

that means in its weak formulation 

f iVul^-^VuVvdx = f fvdx Vf e Wo'P{n \ S). 
Jn Jq 

We introduce the family of diffeomorphisms v^e(x) = x + eX{x) where X is a 
smooth vector field from M.'^ to M^ supported in Q. For e small enough, cp^ maps 
Q into Q. Set S^ = v^e(S) and consider the new state equation in the deformed 
domain 

'- div {\Vu,\f-'Vu,) =f in f] \ S, 

u = on (9fi U S^. 

The corresponding functional is 

J-(E,) = / F(x, u,{x), Vu,{x)) dx + XH\T.,). (9) 

Notice that for e = equation ([8]) reduces to ([7]) and we denote uq simply by 
u. From now on u stands for the unique solution of (jTj) and m^ for ([8]). To 
differentiate ^ we need to show the differentiability of the function e i— )■ n^ at 
zero. If we assume / G L°°{yL) then, by the regularity theory of elliptic equations, 
u and Ue are in ^^^''^{VL \ S) and Wq''^{VL \ S^) respectively. Now let us denote 
by Ue = Us o (fs the transported solution on the fixed domain fi \ S. Since 
Us e W^'^in \ S) (because ip^ is smooth and u^ G W^'P{n \ SJ) the Lemma 4.2 
of [7] gives the differentiability of e i— )■ f/^ from [0, ^o) to Vr^'^(r2 \ S) at zero and 
the differentiability of e i— )■ m^ from [0, Eq) to Wi^^{Q \ S) is obtained via Lemma 
2.1 of |8]. Moreover u' = U' — Vm ■ X in fi \ S, where u' = lim£_^o(w£ — u)/e in the 
distributional sense. For an open set D, taking v,ip G W^'^{D) and tp G C^{D), 
an easy computation of the limit of 

j{Ap{v + t^)-Ap{v),tP), 

as t — )> 0, gives (— div Gy{'Vif),il') where (■, ■) is the duality pairing between T>'{D) 
and C^{D), and Gy is defined by 

G,{Z) = \Vv\P-^Z + {p- 2)\Vv\P-\Vv ■ Z)\/v. 

So the differential of Ap at a point v G W^'''^{D) is given by 

-^iv) = - div(a(Vv.)), y^p G lyi'^(D). 



At u solution of ([7]) the differential of Ap is linear and continuous from W^'^{Q\I1) 
into V'{n \ S) that is 

(9 A (v) 
ov 

The fact that u^ G W^^'^(f2 \ S^) and solves equation ([8]), the differentiability of 
the p-Laplacian operator from W^'P{Q) to V{Q) and the differentiability of the 
transported solution [4 = -u^ o y^^ at e = allow to differentiate the equation ([S]) 
in the distributional sense at e = (see [8J) and we obtain the equation 

-div(G„(VM')) =0 infi\S 

in the distributional sense. This gives an equation satisfied by u' inside fi \ S. It 
remains to find the boundary condition of u'. To this aim, we have to differentiate 
the boundary condition of u^. Due to the particular setting of our problem 
(the boundary operator is the identity and m^ = on dQ U S^), the fact that 
u G iy^'^(r2 \ E) and the differentiability of the transported solution f4 at zero 
imply the differentiability of the boundary condition (see [8j) and we have 

■u' = on dQ and u' = —Vu ■ X on S. 

The fact that u' vanishes on the boundary of Q is due to the compact support of 
the vector field X in Q. So the equation satisfied by u' in the distributional sense 
is 

'-div(G„(VM')) =0 infi\S 

u' = on dn (10) 

u' = —Vu ■ X on S. 

To complete this part, let us check the differentiability at e = of the cost func- 
tion. It is well known (see for example [T])that the length functional '^^(Sg) is 
differentiable at e = 0. The only point to check is concerned with the differ- 
entiability of the map e 1— )■ j^F{x,Ue-,Vui;) dx at e = 0. The smoothness and 
the growth condition on F imply that the map v 1— > -F(-,t', Vf) maps W^'^{D) 
into L^{D) and is differentiable from W^''^{D) into T>'{D) for any open set D. 
Moreover, if u is the solution of equation ([7]), F{-,u, Vu) G W^^'^(i7 \ S) and the 
map e H- F{-,Us,Vus) o fs is differentiable thanks to the hypothesis on u^ and 
U^. Therefore (see Theorem 3.3 of [8]) we have the differentiability of the cost 
function at e = 0. Summarizing, by taking the derivative of the functional ([H]) at 
£ = we get 

— |^^^^(S,) = J (Fy + F, ■ Vu' + div(FX)) dx + X—l^^^n'i^e), 

where u' is the solution of the equation (fTOj) . The derivative ^| _ '^^(Se) that 
appears in the above variation, according to Theorem 7.31 of [Ij, gives 

— \ 



_M\T.,) = f div^ Xdn^ = -{Hs,X). 
Jt, 



As S is countably 7^^-rectifiable, div is the projection of the divergence to the 
approximate tangent hne of S at T-L^-Si.e point of E. 

Unfortunately, the quantity J^{FuU' + Fz ■ Vu')dx is not easily exploitable. 
To overcome this problem we introduce the adjoint state equation 

-div(G„(Vg))=F,-div(F,) in fi \ S 
q = on (9fi U S, 

which has to be understood in the distributional sense 

f (Fu - div{Fz))v dx + f diY{Gu{Vq))vdx = Vf G C,°°(fi \ S). 
Jn Jn 

We are not interested in the regularity of the functions u and q in the whole 
domain Q but only near the optimal set E. The functions u' and q are both in 
W^'^{Q \ S). In the variational formulation of the equation (ITT]) if we take u' as 
a test function, we have 



/ (F„ - div(F^))u' dx+ dw{G^{Vq))u' dx = 0. 
Jn Jn 



(12) 



Let Vf^ and n~ be two sets such that fi = n~^ U ri~ and S C 5^7+ fl dfi~. 
The sets Q~^ and Q~ may be obtained by connecting S to the boundary of Q by 
pieces of smooth curves. The assumption made on S, dQ and / provide sufficient 
regularity for u,u',q so that the Green formula can be applied to f fT2|) . For the 
sequel, we use the following notation: Vu^ stands for the trace on S of Vu 
restricted to ^2"*", ^^ for the trace of the respective normal derivative, -F+ ■ u = 
Fz{x, 0, Vm+) • I'. Similarly Vu~ stands for the trace on S of Vu restricted to Q~ , 
^^ for the trace of the respective normal derivative, F~ ■ v = Fz{x, 0, Vm~) • u. 
Recall also that 

ov 
because m^ = m = on S (i.e the tangential derivative of u^ over E vanishes). 
Let us compute separately the terms of the equation ( IT2l) . Starting by the first 
part, we use the regularity of the boundary of ^2+ to perform the integration by 
parts: 

A+= /" {F^u' - div{F,)u') dx = f {Fuu' + Fz-Vu)dx- f uF^-vdH^, 
Jn+ Jn+ Jdn+ 

where u is the outer normal of Q'^. It is easy to observe that m' = on dQ fl dQ~^ 

and dn+ = s u {dn+ \ {dn u s)) u {dn n dn+) so 

A+ = [ {Fuu' + Fz ■ Vu') dx+ [ ^{F+ ■ u)XudV} - [ u'F^ -udV}. 

Jn+ Jt. c"^ Jdn+\{dnuT,) 



Similarly, taking into account the fact that the outer normal of Q restricted to 
Q'^ n Q^ is — z/, one gets 

A^ = f {FuU + F, ■ Vu') dx- f -^{F; ■ v)Xv dV} + / u'F, ■ v dV} . 

Jn- Jt, (J^ JdQ-\(dQUT.) 



Combining the previous relations and using the fact that the two sets dQ^ \ {dQ U 
S) and dQ- \ {dQ U S) coincide, gives {A = A+ + A') 

A= I {F^u' + F, ■ Vu') dx+ f (^F+ ■ V - ^K ■ '')xiy dV}. 
For the second term, the integration by parts leads to 
B+ = f d\vGu{Vq)udx 

= - f GuiVq) ■ Vu'dx + f u'G^iVq) ■ v dV} 
Jq+ Jdn+ 

where u is the outer normal of Q'^ as in the previous case. It is easily seen that 

B+ = - Guiyq)-Wu'dx- / ^—{Guiyq)^-u)XudH' 
Jn+ JT. (J^ 

+ / u'Gu{Vq)-ydV}. 

Similarly, we have 

B =- [ Guiyq)-Wu'dx+ I ^{Gu{WqY ■u)XudH' 
Jn- is (^^ 

- [ u'Gu{Vq)-vd'H\ 

Therefore, summing up one obtains {B = B^ + B^) 

B = - [ Gu{Vq)-Vudx 
Jn 

+ I (^G„(Vg)- ■ u - ^G„(Vg)+ ■ u'^Xud'H'. 

By the linearity of the function Gu, we get 

/ Gu{Vq)-Vu'dx= / Gu{Vu')-Vqdx 
Jn Jn 

but, by integration by parts allowed by regularity of Q and S, it follows that 

/ GuiVu') ■Vqdx = - f div{Gu{Vu'))q dx + f qGu{Vu') ■ v dV} = 
Jn Jn JanuT. 



because u' is the weak solution of equation fIlOp and q vanishes on 90 US. Finally 
we obtain 



{F^u' + F, ■ Wu') dx = - I {^F+ ■ 



du ^ 



du 



dv 



dv 



■GuiyqY -lAXud-H^. 



To compute the term j^ d\v{FX)dx, we use the divergence theorem thank to the 
reguality of VL. Then 



div(FX) dx 



FXudn 



d-l 







in Jdn 

since X is supported in Q. It follows that 



du~ 



G^{Vq)+ ■u)XiydH^ 



but, by simple computation, we have 

G«(Vg)+ ■ z/ = |Vu+|P-Vg+ ■ z/ + (p - 2)\Vu+\p-\Vu+ ■ Vg+)VM+ ■ z/ 



du^ 



P'^dq- 



dv 

(p-1) 



9z/ 



dv 



- + (p-2) 
p-'^dq^ 



du'^ 



dv 



p-^ / du^ dq^ \du^ 
\ dv dv J dv 



dv 



and also similarly 



Guivq)--v = {p-r 



du 



dv 



p-2Qq- 



dv 



Combining all the computations together we get 
d 



&'-^<^«' 



'\{Hj:,X) 



..|^-=-|^-^ ->--«' 



-{P-K 



du~ 



dv 



p^'^du dq' 



dv dv 



du~^ 



dv 



p-'^du+ dq-^ 



dv dv 



Xvdn\ 



This equality holds for every vector field X, then we derive the following necessary 
condition of optimality: 



X(H._,.) - {^F^ 



9u ^ 

V — 7—Fr ■v] + 



dv 



ip-l] 



du' 



dv 



p-2 Qn Qq- 



dv dv 
10 



du^ 



dv 



p-^du+dq^ 



dv dv 



0. 



We can rewrite this necessary condition of optiniality as 



where the notation (a)^ stands for a^ — a~ . 
We have proved the following result. 



du 



p-2 Qu Qq 
du du 



± 



Theorem 4 Let S he an optimal set in the minimization problem ([2]) and u the 
corresponding solution of the state equation. Assume d = 2, then u satisfies the 
necessary condition of optimality: 



p-2 Qy^ Qq - 



^<-->-(^— (-^)ii:rti)* 



0, 



where v is the unit normal vector o/S, if^ the generalized mean curvature ofH 
and q the solution of the adjoint state equation (ITT]) . 

4 Case of points in M^, d > 1 

In the case of points, some extra difficulties arise because for an equation like 
(ITOl) the gradient and the normal are not defined on points. The equation (ITOl) is 
not the crucial point since we are moving an optimal point xq in the direction of 
the vector field X, the boundary condition of (fTO|) may be writen in the general 
form as u'{xo) = ■§^(^0) if xq is an optimal point. The main difficulty is that 
in the computation we need an integration by part and since the co-dimention 
of the point in M'^((i > 1) is greater than 1, there is a lack of an integration by 
part formulas. The strategy is to study the configurations which are close to the 
optimal one and obtain the optimal configuration as a limit. We consider, for r 
small and positive real number, the set S^ = ip{Br{xQ)) where Br{xQ) is the ball 
centered at the point Xq with radius r and ip is a smooth diffeomorphism from Q 
to Q such that 4j{xq) = Xq. The associated state equation is 



-diY{\Vu\P'^Vu) = f infi\S^ 
u = ondilU K- 



(13) 



For the functional, we consider 

J^(S^) = ^ [ F{x,u{x),Wu{x))dx + ^H'^-\d^r)- 



The factor -^ is in order to avoid the functional to degenerate to the trivial limit 
functional which vanishes everywhere. For r small, 'H'^~^((9Sr) ~ Cr'^~^ where 
C is independent of r and as r — t- the solution of the equation (fT3|) converges 
strongly in WQ^^iVi) to the solution of the same equation defined on i7 \ {a;o}. If 



11 



we denote by u"^ the solution of the equation (I13p for r > and by u the solution 
of same equation for r = (solution on f2 \ {xq}) we may assume pointwise 
convergence of u'^ and Vu^ to u and Vm. Thank to the smoothness of F we have 
also the pointwise convergence of F(-,m^,Vu'') to F[-,u,'Vu) as r — > 0. This 
gives heuristically the scaling factor r^~'^. The same idea for the case of closed 
connected sets in the next section. To simplify the notation in the formulas, we 
will denote the solution of the equation (IT^ by u instead of u^'. Proceeding as 
above, that is transforming the domain by y^e, finding the new state equation and 
the new functional, and taking the derivative of the functional at e = 0, one gets 

— |^^g^((S,,)e) = ^^ / {Fy + F, ■ Vu' + div(FX)) dx - X-^^{Hej:^,X) 

where u' is the solution of the equation 

'-div(Gn(Vu')) =0 infi\S, 

m' = on dVL (14) 

u' = — Vn ■ X on SS^, 

Gu(Vu') is as before and Hq-^^ is the generalized mean curvature of S^. Using 
the fact that xq is optimal and r is small enough (we are in a small neighborhood 
of the optimal point), we obtain ^[^^^^^((Sj.)^) = o(l). 

To overcome the problem of Vu' as in the previous case, we introduce the 
adjoint state equation 

-div(G„(Vg))=F„-div(F,) in fi \ S, 

q = on 9fi U dEr, 

which has to be understood in the distributional sense 

f {FuV-div{F,)v)dx+ [ diY{Gu{Vq))vdx = Vt; G C~(fi \ S,). 

In particular 

/ {Fuu' - div(F^)u') dx+ dw{GuiVq))u' dx = 0. 

By integration by parts, the first term of the equation yields 



/ {Fuu' - dw{F,)u') dx = / {Fuu' + F, ■ Vu') dx - / u' F, ■ Ur dW 



where z/^ is the inward normal of S^. The computation is quite similar to the case 
of closed connected subset of M^ and one gets 



/ 

Jn 



div(G„(Vg))^'rfx = -| ((P-l)|l|^r'|^|^)^^^^'"'- 



f7\Sr JdT,r 

12 



Here SS^ plays the role of E in the two dimensional case. Moreover all the 
quantities vanish in the interior side of E,. then we are interested only on the 
exterior side of S^. We obtain 



n\s, 



{Fuu' + F,-Vu')dx = -f Q^F,-iy-{p-l] 



dv dv dv ) 



Using the above calculation, one can rewrite the derivative of the functional as 



de^ 



X f /du ^ 



du 



dv 



ov ov J 



By a change of variables of type x = iplr, 6), 6 E S"^ ^ we get 

du 






dp 



= 0(1). 



In this notation all the terms of the integrand are evaluated at ilj{r, 9) and J (9) 
is the Jacobian determinant of the function: 6 \-^ iIj{6). It remains to study the 
limit as r tends to 0. We do it in the particular way by letting il){r,6) go to xq 
in a fixed direction as r goes to 0. To express the dependence of the limit on the 
direction ip{6), we use the following notation: v{ip{r, 9)) — )■ i>{iIj{6)) as r — )■ 0; the 
same notation will be also used for other functions in the integrand. This gives: 



[—F,-v-{p 



du 
du 



'■dudq\ 
-—--)XvJde 
dv du/ 



0. 



(16) 



All the terms in the integrand are evaluated at ip{6). The quantity r^~'^{HQj]^,X) 
goes to zero as r goes to zero. In fact if we set fir = r^^'^'H'^~^\-dT,r this measure 
belongs to A4bc and weakly converges to the Dirac mass OJa-i^xg concentrated at 
Xq. Since dimTan(/Xr) = d—1 for all r > it follows that (dimTan(;Ur))/Ur = {d — 
l)fir weakly converges to {d — l)uJd-iSxo. Therefore since dimTan(5o) = < d — 1 
we may apply Theorem [3] with g the constant function d — 1 to have the weak 
convergence of the mean curvature if^^ to the mean curvature Hs^ which is 
identically zero. As a consequence the generalized mean curvature Hgj]^ of S^ 
weakly converges to the generalized mean curvature Hs^ of the point xq. Since 
dH is smooth Hg-^r coincides with the classical mean curvature but to avoid 
confusion with mean curvature of measure in this paper, we keep the terminology 
of general mean curvature. The equality in 0151) holds for every X G C^{Q) and 
every ip diffeomorphism. Again it holds true for X constant in the neighborhood 
of the optimal point and for all ip diffeomorphism satisfying the condition 



/ 



u{ij{9))J{9)d9 = 0. 
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This allows us to write 

, du 



du ^ , 



dv 



p-2 du dq 

---- = const. 



This expression which is evaluated at il){9) is constant for all ip and 9 G S'^^^ . 
This means that it is constant in any direction. Then we have the necessary 
condition of optimality: 



du 



du 



dp 



p-2 du dq 

7—7^ = const, 
dv dv 



Let us consider a particular case of this problem. We assume d = 2 and F = f{x)u 
where u is the solution of the p-Laplacian equation. To express the dependence 
of M on p we denote it by Up instead of u and the same rule for q. Since j9 > 2 we 
want to study the limit as p — )■ 2"*" of the problem. The sequence Up are bounded 
in Hq{Q \ S) then up to extracting a subsequence, it converges weakly to some 
function u. It is easy to see that u coincides with the solution of the classical 
Laplacian that is the solution of equation ([1]) when p = 2. From the adjoint 
state equation I HM we may deduce also that the limit of g^ as p — t- 2+ coincides 
with the solution of equation ([T]) for p = 2. We may then rewrite the necessary 
condition of optimality in the form: 



du 



const. 



dv 
The result proved is summarized below. 

Theorem 5 Let S be an optimal set in the minimization problem ([3]), where 
d > 1 and u is the solution of the corresponding state equation. Then u satisfies 
the necessary condition of optimality: 



du ^ 



du 



dv 



p~2 du dq 

---- = const, 
dv dv 



where v and q are respectively the limit as r ^^ in a given direction of the unit 
normal vector of S,, and the solution of the adjoint state equation (fT5|) . 



The case of points in R is similar to the case of closed connected subsets in 

M. 

5 Case of closed connected subsets in M^ with d > 

2 

Here the strategy is the same. Let S be the optimal configuration. We study the 
configuration which is close to the optimal one and pass to the limit. As in the 
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case of points, we consider a tube H^ = {x & '^'^ '■ d{x, S) < r}. The associated 
state equation is 



- div {\V u\P-^V u) = f infi\S, 
M = on afi U S, 



:i7) 



The procedure is quite similar to the previous case. The corresponding general 
functional is 



^(S, 



1 



H'^-2(5^2) 



A 



^^ F(x,m(x),Vm(x)) + ^^_2^^,_2) 



W-\dT.r 



where 5^ ^ is a (c? — 2)-diniensional sphere of radius r and centered on points of 
E. From the previous computation, we deduce the derivative of the functional: 



d_ 
We 



-L.o-^((s. 



A 



w-^sf-'^) 



{Ha^r^X) 



V-d~\Sfr^) 7aE. V5^ 









where Hqj:,^ is the generalized mean curvature of dT,,.- Remark that all the equa- 
tions are the same as in the case of points in W^. To pass to the limit, we argue 
as in the case of points in W^. First we disintegrate the measure T-C^^^ and get 



^l--^"^' 



A 



■(^as.,X) 



ajd-2( qd- 




E JS? 



du 
dv 



F,-u-{p-l] 



du 



dv 



p-2 



ov ov J 



The measure 'H'^^^{Sf.^^)^^'H'^^'^^Sf^'^ converges weakly to 5^ where x is the cen- 
ter of the sphere S'^"^ . Due to the hypothesis made on data of the problem the 
measure 



(_f,..-(p-i) 



du 



dv 



p-2 du dq 
dv dv 



ld-2i Qd-2\-lnjd-2 qd-2 



weakly converges to the measure 

'-F,-v-{p- 1; 



/du 



\dv 



du 



dv 



vW-^sf: 



p-2 du dq 
dv dv 



v5„. 



for 1H} a.e. z G S. The limit here is computed in a fixed direction 6 G 5"^~^. For 
the curvature part, notice first that the measure /i,. = 'H'^~'^{S'^~'^)~^T-L'^\-dT,r G 
■Mbc weakly converges to the measure "H^lS and (dimTan(/ij.))/ir = dfir weakly 
converges to the measure dl-0'\-T,. The fact that dimTanH^LS = 1 < d allows to 
apply again Theorem |3] to have weak convergence of the mean curvature of yU^ to 
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that of T^^lS and consequently the weak convergence of Hqj]^ to iJ^- Summariz- 
ing all computed results we get 



— I 



^.O^i^e 



\{H^,X) 






(p-i: 



du 



dp 



p-2 



Oh' ov J 



0. 



This relation being true for every vector field X, we get the necessary condition 
of optimality: 



i—F 



{p-l) 



du 



du 



p-2 dudq\ , , ^^ , 



0. 



Theorem 6 Let S he an optimal set in the minimization problem, ([2]), where 
d > 2 and u is the solution of the corresponding state equation. Then u satisfies 
the necessary condition of optimality: 






b-1) 



du 



du 



p-2 Qn Qq 
du du 



X{Hs,u) = 0, 



where u is the unit normal vector of T, in a given direction, iJ^ the generalized 
mean curvature of S and q the limit as r — )■ of the solution of the adjoint state 
equation. 

Remark that this necessary condition of optimality depends on the direction 
9 G S'^~'^ . Those directions are contained in the d — 1 plane which is orthogonal 
to the approximate tangent line to S. 
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